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Introduction 

The  usual  theory  of  shock  wave  propagation  omits 
the  effects  of  viscosity  and  heat  conduction.   It  is  to 
be  expected,  however,  that  these  dlsslpatlve  mechanisms 
will  produce  observable  effects  on  shocks  which  travel 
great  distances.   Since  such  shocks  are  nov;  of  considerable 
importance,  we  have  attempted  to  examine  theoretically  the 
effect  of  viscosity  and  heat  condu.ction  on  them. 

In  this  first  investigation  the  shock  has  been 
assumed  to  be  so  weak  that  all  nonlinear  effects  can  be 
neglected  and  the  usual  acoustic  equations  applied. 
Under  these  assiomptions  the  propagation  of  plane  pulses 
has  been  determined  by  Dr.  J.  Knxidsen,  and  his  results 
constitute  the  first,  and  main,  part  of  this  report.   In 
the  second  part  Knudsen' s  results  are  used  to  find  the 
decay  of  a  spherical  pulse.   This  latter  result  may  apply 
to  spherical  shocks  resulting  from  explosions,  at  great 
distances  from  the  explosion,  when  the  shock  becomes  weak. 

J.  3.  Keller 
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The  Effects  of  Viscosity  and  Heat  Conductivity 
on  the  Transmission  of  Plane  Sound  Pulses 

John  R.  Knudsen 

Fart  I.   General  Discu.ssion 

1.   The  Statement  of  the  Problem 

In  connection  with  the  problem  of  transmission  of 
sound  waves  in  a  compressible  fluid  it  is  of  interest  to 
study  the  effects  of  viscosity  and  heat  conductivity,  since 
these  effects  are  always  present  in  any  real  siti-iation. 
The  purpose  here  is  to  study  these  effects  in  the  case  of  an 
ideal,  polytropic  gas  in  the  semi-infinite  channel  of  Fig.  1. 
V/e  assume,  for  the  purpose  of  simplifying  the  mathematical 
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problem,  that  the  gas  is  initially  at  rest  in  the  channel 
with  constant  pressure,  density  and  temperature  throughout. 

Imagine  that  we  have  constructed  at  the  left  end  of  the 
channel  a  sort  of  fan  or  blovjer  B  directed  to  the  right  into 
the  channel  and  producing  the  same  kind  of  gas  as  is  in  the 
channel.   The  blower  should  operate  in  such  a  vjay  t;?.at  it 
n-enerates  at  any  time  _t  a"  flow  of  gas  vjith  a  prescribed 
velocity  at  x  =  0.   The  introduction  of  this  proscribed  flow 
speed  has  the  desired  effect  of  producing  a  wave  motion  in 
the  channel. 

Assume  the  bloii/er  action  to  cause  a  wave  to  be 
propagated  to  the  right  into  the  undisturbed  gas.   VJhat, 
then,  can  be  said  about  the  form  and  amplitude  of  this 
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wave  and  about  the  general  velocity,  pressure,  density  and 
temperature  distribution  along  the  channel  for  various 
values  of  the  time  t?   For  example,  suppose  that  at  a  certain 
time,  say  t  =  0,  the  flow  speed  at  the  mouth  of  the  channel 
is  suddenly  raised  by  action  of  the  blower  frora  the  value 
zero  to  a  value,  say  1,  which  is  then  maintained  indefinitely. 
A  wave  front  will  then  be  propagated  to  the  right  moving 
into  the  original  lAn^iisturbed  gas  ahead  of  it  and  followed 
by  a  region  of  disturbed  gas  behind  it.   In  the  rather 
trivial  first  approximation  in  which  the  differential 
equations  of  motion  are  linearized  and  in  which  viscosity 
and  heat  conductivity  effects  are  neglected  the  wave  ass^omes 
the  form  of  a  so-called  "square  wave"  moving  x^ith  the  speed 
of  sound  into  the  gas  at  rest  and  followed  by  a  steady  state 
region  x-jith  constant  flow  speed  1.   If  viscosit;.''  and  heat 
conductivity  effects  are  not  neglected  the  form  of  the  wave 
is  distorted  and  the  region  behind  it  is  not  a  steady  state. 
It  ^^fill  be  the  purpose  here  to  investigate  in  a  linear  or 
"acoustic"  approximation  the  effects  of  viscosity  and  heat 
conductivity  on  such  a  wave  motion. 

One  other  simplifying  feature  is  introduced,  i.e., 
the  boundary  layer  effects  along  the  walls  of  the  chann_:l 
will  be  neglected.   Thus  the  problem  is  one-dimensional. 

Of  some  interest  is  a  description  of  the  flow  as  seen 
by  an  observer  located  at  a  fixed  point  x  =  x   for  largo 
values  of  time  t.   This  description  will  show  the  limit  flow 
approached  in  the  channel  behind  tho  wave  front  as  time 
increases  indefinitely  and  the  manner  in  which  this  limit 
is  approached. 

Of  greater   interest   is    the    descriritlon  of   the   wave 

front  as  seen  by  an  observer  moving  with  it.   Accordingly 

WG  introduce  an  observer  moving  along  the  ch,.\nnel  a.t  a  speed 

equal  to  the  sound  speed  of  thu  gas  at  rest,  that  is,  an 

observer  moving  according  to  th:.  law  x  =  x  +  c  t,   T^rom  the 

o   o    "     ■ 

point  of  view  of  such  an  observer  the  velocity  and  pressure 


of  the  gas  in  the  wave  v-jill  depend  uDon  \i   and  X,    the 
coefficients  of  viscosity  and  heat  conductivity,  respectively, 
and  upon  the  time  t.   In  particular,  wo  shall  be  concerned 
with  the  appearance  of  the  velocity  and  pressure  as  a 
functions  of  x   for  large  values  of  t. 

An  even  better  description  of  the  vjave  is  obtained  oj 
an  observer  moving  according  to  the  law  x  =  kc  Jt  +c  t. 
Inasmuch  as  a  fixed  k  interval  will  cover  an  increasingly 
large  portion  of  the  x-axis  such  an  observer  acquires  niore 
information  concerning  the  wave.   This  particular  cbsorver 
is  suggested  by  the  mathematical  analysis  of  the  probloi'ii 
which  is  relatively  simple  in  this  co.go. 

As  far  as  the  nature  of  the  prescribed  V'..;locity  at  the 
left  end  of  the  channel  is  concerned  wo  shall  first  set  up 
the  problem  in  a  general  way  v/ith  an  arbitrary  prescribed 
boundary  condition  at  :c  =  0.   Then  we  shall  discuss  the 
special  case  mentioned  above  in  which,  with  the  gas  initially 
at  rest  in  the  channel,  the  velocity  at  x  =  0  is  suddenly 
increased  at  time  t  =  0  from  zero  to  some  value,  say  1, 
which  remains  constant  for  all  time,   nlso,  as  a  more 
interesting  case  \<!e    take  that  in  vmich  the  velocity  at  x  =  0 
is  increased  from  zero  to  1  at  time  t  =  0  and  then  suddenly 
decreased  to  zero  again  at  some  time  t  =  t, .   We  are 
interested,  in  particular,  in  the  decay  of  the  wave  in  the 
latter  case  as  t  increases  Indefinitely. 

With  respect  to  values  of  u,  and  \  three  situations  will 

be  discussed.   Those  are  the  two  extreme  situations  in  which 

viscosity  and  heat  conductivity  in  turn  do  not  appear  and  an 

intormediate  case  in  which  a  =  ^  =  —  whore  y  is  the  gas 

constant.   (VFhen  necessary  to  give  y  ^-  nuniorical  val.ue  v;e 

assume  the  gas  to  be  air  and  set  y  =  l.li..)   The  cases  [x  =  0 

and  \  =  0,  although  artificial  from  a  physical  standpoint 

present  interesting  mathematical  problems  whoso  conclusions 

at  least  approximate  physical  situations.   The  particular 

value  —  for  ^  as  an  intermediary'-  case  is  chosen  for  two 
Y      A 
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reasons.   For  one  th5.ng  it  serves  to  greatly  simplify  the 
mathematical  analysis  of  the  problem.   Also,  for  y  -  1  •'-!•>  it 

is  quite  close  to  the  ratio  t-  =  -p-. ^  vjhich  Becker  has  shown 

to  be  a  good  approximation  in  a  real  physical  situation. 

2.   General  Results 

In  part  II  will  be  foujid  the  mathematical  solution  of 

the  problem  for  the  cases  a  =  —  ,   a  =  0  and  a  =  oo   v;ith 

Y 
the  two  types  of  waves  descrioed  above  as  viewed  by  the 

several  observers  for  large  values  of  t.   Wo  state  at  this 

point  some  of  the  more  significant  results. 

Let  u,  p  and  p,  all  functions  of  x  and  t,  be  the 

velocity,  pressure  and  density,  respectivelj?-,  of  the  gas  in 

the  channel.   In  the  acoustic  approximation  v;e  let  u  =  6u, 

■o  =  p  +  5p  and  o  =  p  +  5o  and  sot  ud  a  linearized  problem  in 

the  small  quantities  5u,  5p,  and  5p .   The  quantities  p   and 

p   are  the  constant  oressurs  and  densitv  of  the  gas  in  its 

initial  state  of  rest. 

2-1,   The  case  a  =  — 
1 

Consider  first  the  wave  motion  resulting  when  the  flow 
speed  at  x  =  0  is  increased  at  time  t  =  0  from  zero  to  1  and 
maintained  at  that  constant  value  for  all  time  t  >  0.    As 
the  resulting  wave  front  moves  into  the  undisturbed  gas  it 
leaves  behind  it  an  unsteady  state.   The  asymptotic  expression 
for  5u  obtained  by  an  observer  at  x  ~   x  for  largo  values  of 
t  is  given  by 


Actually,  since  5u  is  small  for  all  x  and  t  the  boundary 
value  6u(0,t)  should  be  sm.all.   Because  of  the  linearization 
of  the  problem  Cu  will  at  all  points  of  the  channel  be 
proportional  to  5u'0,t)  and  x^^e  may  consider  the  ratio 
5u{x, t)/5u( 0, t) .   Per  convenience  and  with  no  loss  of 
generality  we.  choose  5u(0,t)  =1.   5u  then  is  essentially 
the  ratio  of  flow  speed  to  the  prescribed  flow  speed  at  the 
entrance  to  the  channel  and  is  dimens ionic ss. 
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(y-1)x  JT- 

( 2 .  01 )       du  ^  1  +  :zr-^T7^  + 


X 


c   is  the  sounci  sveecl   of  the  e;as  at  rest  and  T,  =  ^  is 

proportional  to  \  {and  to  [i   since  [j,  =  — )  and  has  the  dimension 
of  time.   The  second  term  in  the  expansion  is  positive 
indicating  that  behind  the  wave  front  a  flow  speed  greater 
than  that  introduced  at  the  raouth  of  the  channel  is  attained. 
For  increasingly  large  values  of  t  the  flow  speed  decreases 
approachini:';  for  all  observers  x  =  x   the  steady  state  limit 
flow  with  6u  =  1 .   Clearly,  the  smaller  the  value  of  1\  the 

A. 

more  rapidly  is  this  li'nit  flow  approached.   For  T,  =  0  the 

A. 

situation  reduces  to  that  described  earlier  when  viscosity 
and  heat  conductivity  effects  are  neglected. 

An  observer  moving  according  to  the  law  x  =  x   +  c  t 
observes  the  flow  speed  in  the  wave  front.   For  large  values 
of  t  this  speed  is  p;iven  by 


(2.02)   6u(x  +  c  t,t) 
^  ^  o    o    ' 


■  y(y+3)  J  t  *  • 


or  J    in   termr.    of  x   and   t, 


(2.03)       5u,K,t)^l.-i3fi^-Li;£l)(5 
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c^ 
v,  Y+3T  ^  t 


The  significance  of  this  result  is  best  seen  by  considering 

5u  as  a  function  of  x   =  x  -  c  t  for  fixed  large  t  and 

o       o  ^ 
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f  ix(5d  T,  .   The  third  terra,  becuU.sc  of  the  exponential 
function,  is  very  small  for  large  t  compared  with  the  second 
term.   Consider  a  small  x  -interval  centered  about  x^  =  D. 
Fip-.  2  shows  5u  in  this  interval  for  various  values  of  T.  . 


fixed  t 


fixed  T 


\ 


Fig.  2 


For  decreasinr^  values  of  T,  the  x   interval  considered  must  be 
"^  \       o 

decreased  if  (2.02)  is  to  remain  a  good  approximation  for  5u. 

As  1\  decreases  the  curves  become  steeper  approaching  the 

square  wave  form  to  bo  expected  when  viscosity  and  heat 

conductivity  effects  are  not  present.   Also,  for  fixed  values 

of  :c   and  T,  the  velocity  approaches  ^y  as  t  increases 

indefinitely.   This  value  is  halfway  between  the  zero  velocity 

of  the  gas  at  rest  ahead  of  the  wave  and  the  limiting  value  1 

of  the  approximately  steady  state  flow  behind  the  wave. 

In  order  to  stretch  the  x  -interval  and  obtain  an 

o  i* 

expression  for  the  velocity  valid  over  a  greater  range  of 

values  of  X  - c  t  we  introduce  an  observer  moving  according 

to  the  law  x  =  kc  /t'+  c  t.   We  obtain  for  5u  in  this  case 

0     o 

the  asymptotic  expansion 
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(2. Oil-)      ou(kc^Jt'+  c^t^t) 


1  r 


2 


;   l_orf  -^1  + 
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iGysfZ^ 


y(y+3T 


1  -  erf  I 


.r-;^(/T+k)' 


1  e  ■■^"■^ 
71 


J  t 


+ 


or,  again,  In  terms  of  x  and  t 


^    (                   X  -  c  t 
(2.05)   ou(x,t)  -  ^<  1  -cirf  — =:y+     , 


o'  /  .  0-  Sr 


U-  c^t)' 

2T,c^t 
\  o 


j^ 


J£5_ 


[|(Y-1)  J  n. 


YvY^-'' 


.rf  i 


X-C  t| 

X 'I 

"\        C    tf> 

O 


1  c 


2 


T 
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In  a  small  k  Interval  c:;ntered  about  k  =  0  we  consider 
6u  for  large  fixed  t.   Pig.  3  shows  Su  for  various  values 
of  T  .   For  decreasina;  values  of  T.  the  curves  become  steeper 


Fig.  3 


approaching  the  square  wave  limit  situation.   For 
indefinitely  largo  t  the  velocity  5u  becomes 


9. 


1-  erf 


/^^ 


v/hlch  in  any  k  Interval  depends  on  T^ .   Thus  more  information 
is  obtained  by  our  present  observer  than  by  the  earlier  one 
for  which  the  ultimate  value  of  5u  turned  out  to  be  sj.mply  p* 
over  any  snail  x  -Interval. 

In  the  event  that  after  a  relatively  small  finite  time 
t,  the  velocity  at  x  =  0  is  suddenly  reduced  to  zero  a  second 
wave  front  is  propagated  to  the  right  in  the  channel.   It 
moves  into  the  wake  of  the  first  wave  front  causing  it  to 
decay  or  die  out.   From  the  better  vlevjpoint  of  the  observer 
X  =  kc  /t"+  c  t  we  obtain  for  Su  the  asymptotic  expansion 


( 2 . 06 )      5u 


/27iT.t 

A 


I,' 


kt. 


2T 


/2^ 


y(y+3T  J7t 


1  -  e 


Yt. 


\ 


-T7^(Jt  +  k)2 


^■\ 


1  -  erf 


k 


x  -  c   t 
or,    with  k  =  


(2.07)       5u(x,t) 


t^  x  e 


2c^T, t 
o   \ 


C     tf7Ti'L\  t 
o  \ 


Yt. 


Y(  Y  +  3)      VTT    \ 


m 


\ 


J2L 


J+cVt    IT-    ,  . 

oX  XJt       ,^-f>l 


—  x-c^t 

j: 2_ 

-A      c^/t 


+  ... 


For   a   fixed  k   interval   the    first   term  predominates  and 
-1/2 


6u 


0  as    t 
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The  asymptotic  expansion  for  the  pressure  in  this  same 
instance  is  given  by 


(2.00)      p   =   p^+  6p 


-^o 


1  + 


'"2Y^-3Y+h  ,  2y^+3y-^     X    1  V 


(X-  c^t)' 

2c^T-t 
o   \ 


Ti 


!-Y 


ir^ 


c^ti 


T^wrr 


+ . . . 


-]/2 

v/hence    6p  — =»•  0  as    t    ''       also. 

Of  interest  too  is  the  flow  behind  the  second  wave 
front.   For  an  observer  at  fixed  x  =  x   the  velocity  is 

o  "^ 

given  by 

(2.09)      5u(x,t)  ^ 


M-CqJy^  t- 


Thus  the  gas  has  a  negative  velocity  in  the  v;ake  of  the 
second  wave  front.   The  limitiiig  state  of  rest  is  approached 

as  t  -^^  .   This  contrasts  vjith  the  approach  to  the  steady 
state  with  5u  =  1  behind  the  first  wave  front  when  the 
second  wave  is  not  present.   The  approach  in  this  case  is 

-■3/2 

proportional  to  t  -^^  , 


2-2.   The  case  a  =  0 

In  this  somewhat  artificial  case  in  which  viscosity  is 
absent  but  heat  conductivity  present  the  results  are  similar 
to  those  of  the  preceding  case.   Ineod  in  determining  the 
flow  behind  the  second  wave  for  fixed  x  and  large  t  the 
result  is  Identically  that  given  by  (2.09). 

Passing  over  the  observer  x  =  x  +c  t  and  considering 

only  X  =  kc  >rt  +  c  t  we  observe  the  tvso   X'jave  combination  to 
have  flow  speed 
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(2.10)       5u  ^ 


2T^  (5Y-l)kt, 


X 


1        c 


"\ 


f^^fvTt 


8i  2TtY(  Y-TTT 


VfrT 


-.^{ 


-[.[.3f-^  .99  #-.75|^; 


1  -erf   1.02 


^J 


,  •<  sm  ! 


sin  {  .03  ~+  .i-^5  ^+  .73  ~+  .20 


.56 


-   e 


sm 


t    .     ,  ^  k-Tt  .     ^.   k 


v2 


-] 


.Ii5  Fif-+  .ij-5  ~+  .73  ^+  .20 
■^x  X  "X  ^ 


Compare  this  with  (2,06). 


2-3.   The  case  a  =  oo 

In  this   other   extreme    situation  of  i'leat    conductivity  by- 
no  viscosity   the   matheraatics    simplifies   sor:iovrhat,    the 
differential   equations   being   of   lower    order.      Again   the 

results   are  not  linlike    those    of   the    case    a  =  —  . 

Y 
For  large    t   the    observer  x   =  x    +  c   t   finds   for   the 

single   pulse  wave   front    the   flow    speed 


(2.11) 


Cu 


11 


J^ 


3 

7J 


X 


O     [I 


T 
t 


(T   =  — ■^^)  which  compares  with  (2.03). 

p  c 
•^o  o 

For  the  observer  x  =  kc  Jt  + c  t  the  situation  is 

o     o 

described  by 
(2.12) 


5u  -v^  i  "fl  -  erf   -^ 


Jtl 
^    "  2T   [  T 


u. 


8^2^" 


+  . .  . 


which  compares  with  (2.01l ). 
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Part  II.   Mathematical  Analysis 

3.   The  Differential  FJquations 

The  differential  equations  for  a  one -dimensional 
compressible  gas  flow  with  viscosity  and  heat  conductivity 
terms  present  are 


fp^    +  (pu)^  =  0 


(3.01)  ; 


(pu),  +  (pu   +  p  -  l^u  )   =  0 


XX 


^  [p(|u  +e)  ]^  +  [pu(|u  +i)  -  M-uu^-  Xc^T^]^  =  0 


where  p(x,t)  is  density,  u(x,t)  velocity,  p(x,t)  pressure, 
e(x,t)  energy,  i(x,t)  enthalpy, .T(x,t)  absolute  temperature 
and  c   specific  heat  at  constant  volume. 

We  define,  for  convenience,  the  quantity  9, 
proportional  to  T,  bj^  p  =  pO.   Empirically,  we  have 


(3.02)    e 


9 


Y 


-  1 


1  =  —i. 


'9 


@ 


Y  -1   ' 


c^(Y-l) 


and  equations  ( 3 • 01  become 


(3.03) 


rp^   +    ipu)      =   0 


(pu)^  +  (pu  +  p9  -  M-u^)^  =  0 


^p^h^^^'h^  'p''^h^-'T^^-^''''yr-^^x  = 


0 


¥e    introduce  an  acoustic  approximation  by  setting 


(3.  Oil-)     u  =  6u   , 


0  =  ^0+  5p   ,    9  -  9^+  0© 


where  o   and  9  are  the  constant  density  and  temperature  of 
"^o      o  ^^         - 

the  gas  at  rest.   Making  these  substitutions  in  (3.O3)  and 
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neglecting  non-linear  term?  we  obtain  differential  equations 
for  the  small  quantities  6p(:<,t),  5u(x,t)  and  5©(x,t),  i.e., 


(5p)^  +  Po(6u)^  =  0 


(3.05) 


Pq(5u)^  +  Po(^^®)x  ""  ®o^^P^x  "  l^(S^^ 


XX 


p^loe)^  +  (y-1)Po©o(5u)^  =  X(50) 


XX 


[j-.   Initial  and  Boundary  Conditions 

Th3  differential  equations  (3.05)  are  to  bo  solved 

subject  to  certain  initial  and  boundary'"  conditions.   For 

simplicitjr  we  assume  the  initial  state  to  be  that  of  the 

original  undisturbed  gas  v/ith  constant  density  and  temperature 

p   and  9  .   Thus 
•^o      o 


(l+.Ol) 


5u(x,C)  =  5p(x,0)  =  6G(x,0)  =  C 


At  the  left  ena.   of  the  cha:miol  vjc  impose  two  boundary 


conditionc 


one  on  6v.   and  one  on  6Q,      As  it  t'orns  out  these 


are  sufficient  for  the  evaluation  of  constants  of  integration 
occurrino;  in  the  solution  of  the  differential  equations. 
The  absence  of  op,   from  the  differential  equations  seems  to 
to  indicate  that  no   condition  should  be  imposed  on  5p  at  x=C, 

As  the  boundary  condition  on  6u  vje  have  the  proscribed 
velocity  at  the  mouth  of  tho  chann.el  due  to  the  controlled 
action  of  the  blower. 


(lj-.02) 


5u(o,t)  =  Tj(t; 


As  a  physically  reasonable  condition  on  the  temperature  we 
insist  that  there  bo  no  transfer  of  heat  into  the  entrance 
of  the  channel,  i.e., 


(i|.03) 


6fe,^(0,t)  =  0 


Ik- 


As  far  as  the  right  extreraity  of  the  channel,  .:c  =  oo , 
is  concerned  the  natural  conditions  are 


(^.01+) 


5u(  00  ,t;  =  5p(  00  ,  t  ) 


5G(co  ,t)  =  0 


The  latter  conditions,  however,  are  taken  care  of 
automatically  by  the  method  of  solution  of  the  differential 
equations . 

5»      Solution   of   the    Initial-Boundary  Val^ie^Problorti 

The  solution  of  equations  (3«05)  with  initial 
conditions  (k.Ol)  and  boundary  conditions  (;4..02),  (k,03) 
and  (ij-.OLi-)  is  accomplished  bj   means  of  a  Laplace  trans- 
formation consisting  of  multiplying  each  term  of  equations 

-s  t 
(3»Cp)  by  e    and  integrating  between  the  limits  t  =  0 

and  t  =  CO  .   The  functions  6u(x,t),  5p(x,t)  and  50(x,t) 

are  transformed  into  w(x,s),  r(x,s)  and  h{x,s)  respectively 

Equations  (3»05)  with  the  aid  of  the  initial  conditions 

([). .01)  are  transformed  into 


sr  +  p^w^-  0 


(5.01) 


J  p  sw  +  p  h  +  ©  r   =  uw 
-^  '  o     'ox    ox    •  XX 


p  sh  +  (y-I'^o  6  w  =  \h 

'  O       ^  '     ■OCX       XX 


and  the  boundary  conditions  (Ji.02)  and  {U.,0^)    into 


(5.02) 

and 

(5.03) 


w(0,s)  =  U(  s) 


hjO,s)  =  0 


The  transformed  boundary  value  problem  defined  by 
(5.01) J  (5.02)  and  (5.03)  can  be  expressed  in  terms  of  the 
single  unknown  function  w(x,s)  by  eliminating  r  and  h.   We 
obtain 


15. 


2  3 


(5.0L)   ^(Po©o+l^s)Wxxxx-  Po^WPoQo^!^^^^^)^^xx^Po^^^'  =  ^ 


as  the  differential  equation  for  w(x,s)  and 


(5.05) 


w(  0,  s)  =  U(  s) 


and 

(5.06) 


P^s   U(s) 
w   (C,s)  =  ^ 


as  the  boundary  conditions. 

The  appearance  of  (5. Oil-)  is  simplified  if  we  write  it 
in  the  form 


(5.07) 

xvhere 

(5.03) 


^xxxx-  (^^-^^  °^'^\x^  '^^'^  =   0 


a(s) 


b(s) 


■S^J-H-p-Tr^i  YPoV'^^-^'-^-J  ^^PoV!^^-'^^)  -l|.Xs(p^G  +us)| 
■  o  O  '     ^  -^ 


in  which  the  +  sign  refers  to  a(s)  and  the  -  si.3;n  to  b(s) 
The  solution  to  (5.0?)  satisfying  conditions  (5.05) 
and  (5.06)  is 

(5.09)   w(x,s) 


r    Po- 


TJ(s) 


PoVl^^^ 


b- 


2     PQ^ 


a  - 


2   ,  2 
a  -  0 


e-a(s)x  ^ 


PoV*^^   .b(s): 


2  ,  2 
I  -  0 


or 
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(>.10)        w{x,s) 


n(s) 

2 


YP,^©j^+  us  -Xs 


1    -a' 


s)x 


J  (  YPo^o"^  '^^  "^ '^^^     -  l4-Xs(p^e^  + tis)  J 


1   +  — 


YP    ©,.,  +  US  -  \s 


J(yPo&o+  t^s+  Xs)---  ]+\s(p^©^+  lis)    . 


-bC  s)x 


J 


(^.11) 


Tho   notation   is    soraewhat    simplified   if  wc    introduce 

T 


■'00 


■p.        YP„^)^    '      ^        'XS    *         ^        T^        X         • 


The  quantities  T^^  and  T   are  proportional  to  X  and  [x, 
respectivolj'-,  and  have  the  dimension  of  time.  d"   and  a  are 
dimansionless.   The  solution  (5.10)  becones 


(5.12)        w(x,  ^) 

■'x 


■"x 


1  - 


1+  (a-DO" 


.    -a(Tf  )x 


waer 


T 


X 


/ri+  (a+Dcr]"^  -  -cr(i+  YO'J'i 


+ 


1  + 


1+  (c-Do- 


e         ^ 


'-     /Ti  +  ( a+i)  cT]  ^  -  ^(  1  +  Ytt  or)  ^ 


(5.13) 


c    T- 
o    X 


^n 


i^^~aT  (^^  (ci+Dcrt  J[i+  (a+i)(r]2-  ^cr(i+  ycct)}   . 
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(c   =  Jy^o  ^^   ^^^    sound  speed  of  the  gas  of  rest.) 

The  next  step  is  to  obtain  the  desired  velocity  function 
5u(x,t)  from  its  Laplace  transform  v;(x,  7^).   This  is 

accomplished  by  means  of  the  inversion  integral 


(5.lL^)      5u(x,t)  =  tIi  /  ^^(^^'f-) 


T~^ 

T 
Br       ^        ^ 


where  the  path  of  integration  is  the  usual  Bromwich  line. 

It  is  clear  from  (5*12)  that  the  number  and  location  of 
branch  points  and  poles  of  the  integrand  (5«l2)  depend  upon 
the  value  of  c.  and  \ipon  the  prescribed  function  whose  trans- 
form is  U(^).   Thus,  in  order  to  continue  to  an  explicit 

solution  we  must  specify  a  and  U(t).  As  already  mentioned 
we  V7ill  discuss  three  values  of  a,  i.e.,  —  ,  0  and  00.  V/e 
have  mentioned  also  the  two  cases  we  shall  consider  for  'J(t). 

6.   The  Inverse  Transform 

The  actual  evaluation  of  (5.1[|.)  for  the  several  cases 
mentioned  is  a  lengthy  process.  We  3;ive  here  only  the  more 
significant  steps. 

For  the  case  a  =  —  equation  (5 •12)  becomes 

(6.01)  w(x,^)  =  arTj(^).  ^-1-^e  °  ^ 


C_T 
+  ■::rr — —, — ^-,  x  ,  -   e 


cr  [  -y  +  (  Y  - 1 )  (S" 


o-\  /I+a' 


(Wo  see  here  the  reason  for  the  choice  of  —  as  a  value 

Y 
for  a.   It  is  just  the  value  of  a  tdat  makes  a  perfect  square 

of  the  quantity  under  the  radical  signs  in  (>.l2).   Thus  it 

eliminates  possible  branch  points.) 


18. 


77/  CT 


Considering  IT(^)    wo   havo    tv;o    situations.      If  U(t)    =  0 
for   t  <    0  and  U(t)    =   1   for   t   >   0  then  U(^)    =  7^  .       If 
IT(t)    =   0  for   t   <    0  and   t   >   t^^  and  y(t)    =  1   for   0  <   t  <   t-j^ 

T     /■  -A 

then  U(?p— )    =  "Tp  I  1  -e  j    .      In  both   cases    it    contributes 

a  pole    to   the    integrand   of    (5.12). 

T 
With  U(~)    =  ~      (6.01)    may  be  written  as    the    sirm  of 
-\  " 

two   integrals 

t   ^  X        1 

r         if       '^    1        T~^-T^-^"J^ 

(6.02)      i^(x,^)    =   2-TtT  J     T-rtY~W 


^        '°^^  dcr 


Br 


and 


t    _          X 
rr—  cr  -  -;: 


/^  C      T  r— — — 

Y  "X  o^X  /l+?   ,^ 

^  Br 

These  are  to  be  evaluated  for  each  of  three  observers  x  =  x  , 

X  =  X  +  c  t  and  x  =  kc  JT + c  t.   Wo  cannot  obtain  explicit 
00  00  ^ 

values  for  these  integrals  for  all  values  of  t.   Hovjever,  we 

can  find  asymptotic  expansions  v;hich  are  valid  for  large 

values  of  t  and  this  is  ^^rhat  we  are  mainly  interested  in. 

For  the  observer  x  =  x   the  expansions  are  easily 

obtained  by  choice  of  a  proper  path  of  integration  equivalent 

to  the  3r omw i ch  line. 

For  X  =  x^ +  c^t  the  integral  I,  assumes  the  form 
00  "^1 


I  =^.Jc^(cT)e*^(^)d.(r 


v/hich  by  the  saddle  point  method  of  steepest  descents  has  an 
asymptotic  expansion 
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1^2^  sj  I  f"  ( <r^  )TJt 


where   CT  is  the  saddle  point  and  p  the  direction  of  the 

path  of  integration  through  it. 

The  integral  I„  does  not  yield  to  the  method  of 

steepest  descents  since  the  saddle  point  (f     of  the  real 

part  of  f((7')   is  a  pole  of  the  integrand.   The  expansion 

of  Ip  is  obtained  by  an  appropriate  deformation  of  the 

Bromwich  path  and  reduction  to  standard  integrals. 

For  the  observer  x  =  kc  /t  + c  t  the  integral  I, 

o     o  '^     1 

assumes  the  form 

I  =^  J'^(a')e*^'^^^^g(^)  do-   . 

3r 


The  asymptotic  expansion  in  this  case  requires  a  modified 
saddle  point  method  based  again  on  the  saddle  point  of 
the  real  part  of  f(<J')  .   The  expansion  one    obtains  is 

(j)  ( or    )  e  "       G   ■" 

-r  S 


iJz^,  7TF~((rjl7T 


!g'(GrJ 

1  -  erf  ^ 


2if''(  a;)l 


(For  brevity  the  development  of  this  is  omitted.) 

The  integral  Ip  again  resists  the  saddle  point  method 
but  is  handled  by  a  proper  deformation  of  the  3romwlch  path, 


20. 

In  the  case  TJ(^)  =  -^  M  -  e   '^   j  the  integrals  are 

not  essentially  different  than  in  the  preceding  cases. 

In  th'3  situation  closcribod  "oj   a   =   Q   the  integrand  of 

('^.li;)  is  somewhat  more  complicated  than  in  the  case  o,  =  — 

Y 
because  of  presence  of  more  branch  points  and  poles. 

Hov-fcver,  the  saddle  point  methods  described  above  suffice 
to  evaluate  the  integrals. 

V/hon  a  =  CO   something  different  occurs.   The 
differential  equation  for  w(x,s)  reduces  to  second  order 
and  one  boundary  condition  has  to  be  dropped.   Rather 
naturally  we  hold  on  to  the  one  prescribing  the  velocity 
of  the  blower  and  drop  the  one  governing  the  heat  transfer 
at  X  =  0. 

The  one  integral  arising  in  this  case  is,  for  all 
observers,  evaluated  by  proper  deformation  of  the  Brom.wich 
path,  and  reduction  to  standard  integrals. 
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